Polarization properties and dispersion relations for spiral resonances of a dielectric rod 
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Dielectric microcavities based on cylindrical and deformed cylindrical shapes 
have been employed as resonators for microlasers. Such systems support spi- 
ral resonances with finite momentum along the cylinder axis. For such modes 
the boundary conditions do not separate and simple TM and TE polariza- 
tion states do not exist. We formulate a theory for the dispersion relations 
and polarization properties of such resonances for an infinite dielectric rod of 
arbitrary cross-section and then solve for these quantities for the case of a 
circular cross-section (cylinder). Useful analytic formulas are obtained using 
the eikonal (Einstein-Brillouin-Kcllcr) method which are shown to be excellent 
approximations to the exact results from the wave equation. The major finding 
is that the polarization of the radiation emitted into the far-field is linear up 
to a polarization critical angle (PCA) at which it changes to elliptical. The 
PCA always lies between the Brewster and total-internal-reflection angles for 
the dielectric, as is shown by an analysis based on the Jones matrices of the 
spiraling rays. 
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1. Introduction 

There has been a great deal of recent interest in cylin- 
drical and deformed cylindrical dielectric resonators for 
micro-laser applicationsiii2i2i4£ From the theory side 
there is a particular interest in the deformed case, as in 
this case such resonators are wave-chaotic systems and 
can be analyzed with methods from non-linear dynamics 
and semi-classical quantum theory. Analysis of the reso- 
nances and emission patterns from such systems has fo- 
cused exclusively on the scalar Helmholtz equation which 
describes the axial component of the electric (TM mode) 
or magnetic (TE mode) fields for the case of resonant 
modes with zero momentum in the axial direction (z- 
direction). For this case (k z = 0) the polarization state 
is unchanged by boundary scattering and the non-trivial 
ray dynamics in the transverse plane does not affect the 
polarization state of the resonant solutions. However it is 
interesting to consider the solutions of the wave equation 
for both cylindrical and deformed cylindrical dielectric 
rods with k z =^ 0, since in this case the boundary scat- 
tering couples the electric and magnetic fields and there 
no longer exist TM or TE solutions with a fixed direction 
in space for one the fields. We refer to these non-zero k z 
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modes as "spiral modes" and note that elastic scattering 
from such spiral resonance modes has been measured pre- 
viously by Poon et al. Ref. 0, the polarization properties 
were unfortunately fully explored in those experiments. 

The authors did however measure a systematic "blue- 
shift" of the resonance mod es with tilt angle, which has 
been predicted by Ref. I7l8l and which we will derive be- 
low. The modes we study here are similar to the hy- 
brid modes known in the study of optical fibers, where it 
is also well-known that there are no simple TE, TM or 
TEM-like modes, but more complex vector solutions are 
necessary. Our emphasis differs however in several ways: 
1) we are interested in uniform dielectric rods, not the 
variable index profiles typical of optical fibers. 2) We are 
interested in modes which are not totally internally re- 
flected so we can study the nature of the polarization of 
the emitted radiation in the far-field. 3) We are primar- 
ily interested in resonances of uniform rods with cross- 
sections in the range of tens to hundreds of /zm, so they 
are strongly multi-mode and can be treated within the 
eikonal (semi-classical) approximation. 

Here, we will model finite resonators as infinite dielec- 
tric rods, neglecting the end effects in the z-direction; 
we therefore formulate the vector wave equation and the 
necessary boundary conditions for an infinite rod of arbi- 
trary cross-section. Various approximations are possible 
to treat end effects when they are relevant, but we will 
not explore them here. We define the resonant solutions 
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Fig. 1. Schematics of a spiraling ray in a rod with ar- 
bitrary cross-section. The ray can refract out and the 
polarization can be defined in the far-field with respect 
to a local coordinate system. For a general deformation, 
the polarization will change with the azimuthal angle <f>. 
Here, we will focus only on the circular cross-section, 
where the polarization is constant (with respect to cf>) 
when referred to the local coordinate system oriented 
along and perpendicular to a ray in the far-field. 



(quasi-bound modes) of such a system and write down a 
general formalism which can be used to obtain exact nu- 
merical solutions for the vector resonances. We also show 
how the corresponding solutions can be used to derive the 
spatially- varying polarization state of the emitted radia- 
tion in the far-field. We then study in detail these equa- 
tions in the case of a circular cross-section (cylinder) for 
which a great deal of analytic progress and physical in- 
sight may be obtained. The analysis of the non-circular 
case for which wave-chaotic polarization states are pos- 
sible will be published elsewhere^ 

Our goal is to relate the polarization state in the far- 
field to the projected two-dimensional ray motion in the 
plane transverse to the z-axis. (See Fig.^l One can have 
spiral resonances which range from motion along the di- 
ameter of the rod in the transverse plane (bouncing ball 
type modes) to whispering gallery modes which circulate 
around the perimeter of the rod as they spiral along it. 
We will discuss how the polarization properties of the 
resonances vary as the angle of incidence in the plane 
(sin%) and the spiral (tilt) angle (6) with respect to the 
x — y plane varies. (In Fig. [5] A) we introduce the rel- 
evant geometric parameters.) It should be noted that 
due to the curvature of the boundary, even modes which 
are totally-internally reflected according to geometric op- 
tics do emit by evanescent radiation into the far-field and 
their polarization fields can be obtained from exact so- 
lution of the wave equation, although experimentally it 
may be impractical to measure their weak emission far 



above the critical angle. In the simplest case of k z = 
{9 = 0) one finds pure linear polarization in the far-field; 
in addition, the resonant energies are just those of the 
two-dimensional problem of a dielectric disk. In this two- 
dimensional case the resonant energies in the semiclassi- 
cal limit are determined by the optical path length as 
well as by the phase shifts due to caustics and reflections 
at the boundary^ These boundary terms in the semi- 
classical limit correspond simply to the TIR phase shifts 
for TM and TE scattering off a plane dielectric interface 
when sin x > sin Xc = 1/n (here n is the index of refrac- 
tion of the rod surrounded by air); if sinx < sinx c = 1/n 
there is zero phase shift but just a loss (imaginary part 
of k) given by the Fresnel scattering coefficients. For the 
spiral modes (k z ^ 0), the boundary terms have a new 
character corresponding neither to the TE nor to the TM 
Fresnel scattering and new phenomena can occur, such 
as a non-zero phase shift for modes which are not TIR. 
We derive below the generalization of the Fresnel scat- 
tering coefficients for spiral modes of the cylinder using 
the vector eikonal method. We find that the angle at 
which a non-zero phase shift sets in is always between 
the critical angle and the Brewster angle and coincides 
precisely with the onset of elliptical polarization in the 
far-field. We call this new key quantity the polarization 
critical angle (PCA) 

In Section of the paper we set up the relevant form 
of the vector wave equations for the infinite dielectric 
rod and formulate the boundary conditions for the quasi- 
bound modes (resonances). In Section we discuss how 
to extract the far-field polarization of these quasi-bound 
modes. All results are for the general case of arbitrary 
cross-section of the rod. In Section0]we specialize to the 
case of the dielectric cylinder and reduce the resonance 
problem to a simple root-finding problem. This equation 
is exact and is shown to yield the systematic blue-shift in 
Rcf. 6. In Section[S]we reformulate the same problem us- 
ing the eikonal (EBK) method which yields a simpler an- 
alytic formula for the resonances and allows a statement 
of the polarization problem in terms of generalized Fres- 
nel coefficients. The exact and EBK resonance wavevec- 
tors are shown to agree quite well, down to small k. In 
Section|S|we restate the polarization problem in terms of 
Jones matrices and thus derive the internal polarization 
state and the farfield polarization of the resonances in the 
semi-classical limit. Both the Jones and EBK formula- 
tions are shown to yield the same answer for this quantity 
and to agree with the exact results to a good approxima- 
tion. Finally, the origin of the polarization critical angle 
is explained. 



2. Wave equation and resonances for the infinite 
dielectric rod 

For electromagnetic fields in free space interacting with 
uniform dielectrics, Maxwell equations reduce to the vec- 
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Fig. 2. A) Coordinates used for ray dynamics in a rod with arbitrary cross-section. The ray can refract out and the 
polarization can be defined in the far-field with respect to a local cartesian coordinate system tilted so that one of the 
axes is along the propagation direction. We define rj as the angle of incidence in the plane of incidence; \i t he angle of 
incidence projected into the plane; cr, the projected far-field angle; the tilt angle measured from the cross-sectional 
x — y plane, given by tan# = k z /ji; a, the external tilt angle tana = k z /-f2', B), C) Schematics highlighting the 
trigonometric relations among the inside and outside quantities. Note that a can be found through the application 
of Snell's law, rising — sin a. 



tor Helmholtz equation 



E(x,y,z) 
B(x,y,z) 



0, 



(1) 



where n is the uniform index of the rod and n = 1 out- 
side the rod. Thus n differs from unity in an arbitrary 
closed simply-connected domain 3D in the x — y plane 
for all values of z (the rod is infinite in the z-direction). 
The translational symmetry along the z-axis (see Fig.^) 
allows us to express the ^-variation of the fields as 



E(x) =E(x,y)e 



ik-2 



B(x)=B(x,y)e 



-ik z 



(2) 



and henceforth the vectors E, B will refer to the x,y 
dependent vector fields just defined. With this ansatz, 
we can show that the most general solution of the six- 
component vector wave equation for this problem is de- 
termined by E z and B z components alone; the perpen- 
dicular fields are given by linear combinations of these 
two scalar fields and their derivatives ~ Hence we must 
solve the two-component scalar wave equation 



with 7 2 = n(xfk 2 - k 2 z 



(3) 



where we have introduced the reduced wavevector 7 which 
is the wavevector associated with the transverse mo- 
tion. The complication of solving this remaining two- 
component Helmholtz equation stems from the fact that 
the two fields E z , B z are coupled through the boundary 
conditions. Four independent boundary conditions are 
found through the application of the general Maxwell 
boundary conditions: 



E zl = E z2 or 8 t E zl = d t E z2 (4) 

B zl = B z2 or d t B zl = d t B z2 (5) 

-^d n B z i yd n B z2 = — I — o 1 I d t E z i (6) 

7i 72 \7i 72/ 

—^dnEzi %~d n E Z 2 = + ( — I 1 ) d t B z i. (7) 



7i 



I2 



7f 75 



Here, 1 and 2 refer to inside and outside solutions. We re- 
cover the familiar special case of two-dimensional modes 
when we take k z = 0; in this case the boundary condi- 
tions can be satisfied with either B z — (TM solutions) 
or E z — (TE solutions) i.e. we get spatially uniform 
eigenpolarization directions. For the TM case we find 
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Ei 



and 



B 



—zn k 



-d v E, 
d x E z 



and for the TE case we find B i = and 



7 



, , —d y B z 
2 *l d x B z 



(8) 



(9) 



In both cases the equations Eqs. ©-© completely de- 
couple and we have only to solve the scalar Helmholtz 
equation for E z or B z with the boundary conditions of 
the continuity of the field and its normal derivative on 
the boundary. In both cases the electromagnetic field is 
linearly polarized in the far-field, either in the direction 
parallel to the rod axis (TM) or perpendicular to it (TE). 
We will now focus on the k z ^= modes which we will re- 
fer to as spiral modes. We will not consider the extreme 
case where k = k z (TEM modes). 

Our interest is in the dielectric rod as a resonator, i.e. 
as a device for trapping light. Experiments on resonators 
fall into two broad categories, and the presence of quasi- 
bound modes are manifested differently in these two sit- 
uations. One can measure elastic scattering of incident 
laser light from such a rod, as was done in Ref.^and focus 
on the specific wavevectors at which one observes scat- 
tering resonances. For this case the linear wave equation 
that we are studying provides an exact description. One 
can also imagine the rod containing a gain medium and 
when pumped emitting laser light into these spiral reso- 
nances^ In such a case, the linear wave equation is not 
an exact description; however for high-Q resonances typ- 
ically the resonances of the passive and active cavity are 
very similar^ In order for the laser to emit specifically 
into spiral modes, there would have to be some mecha- 
nism to suppress lasing of planar (k z = 0) modes; one 
could imagine doing this with a small seed pump which 
is tuned to the frequency of a spiral mode and pushes 
it above threshold before all the other modes. It is pos- 
sible to describe these two different physical situations 
(elastic scattering and lasing) using appropriate bound- 
ary conditions on the linear vector Helmholtz equation. 
To be precise, we assume that the rod is bounded by the 
interface dD given by (see Fig. 

dD = R(z,<j)) V z£M, 4>€[0,2ir}. (10) 

We are typically interested in boundaries which are 
smooth and not too far from a circle, hence it is natural 
to expand the internal and external solutions E z and B z 
of Eq. © in cylindrical harmonics^ 



E< 
B< 



E> 
B> 



E 

1 — — o* 

oo 

E 



H-(lir) 



im<f> 



(ii) 



) H+{l2T) + fr) H m { l2 r) 



vm<j> 



(12) 



where are the Hankel- functions, H representing an 
incoming wave from infinity and H + an outgoing wave. 

To describe the scattering experiment we simply ap- 
ply the boundary conditions Eqs. which will 
connect these interior and exterior solutions by a set 
of linear equations for the coefficients a m , £ m , . . . , d m . 
These linear equations will have solutions for any inci- 
dent wavevector k and will define a scattering matrix 
for the fields which can be used to calculate the intensity 
scattered at any given far-field angle. An example of such 
a calculation for the cylinder is given in Fig. EI the values 
of kR (R is the cylinder radius) at which rapid variation 
is observed are the resonant wavevectors for which the 
incident light is trapped for long periods. However the 
precise pattern of radiation in the far-field in this case is 
determined by both the scattered and incident radiation 
and is not representative of a lasing mode for which there 
is no incident radiation. To determine the resonances cor- 
responding to emission from a source it is conventional to 
use the Sommerfeld or radiation boundary conditions; in 
this case we would set the incident waves from outside to 
zero (the coefficients 6 m , $m in Eq. iJH}) an d still impose 
the boundary conditions of Eqs. Q-Q. The resulting 
linear equations for the remaining coefficients would not 
have any solutions for real k (since current is not con- 
served), but would have discrete solutions for complex 
k values; those solutions are known as the quasi-bound 
modes of the problem. It can be shown that the discrete 
complex solutions of this problem correspond to the poles 
of the S'-matrix of the current-conserving problem when 
that S- matrix is continued to complex values of k (see 
e.g. Ref. Il0ll3^) . Re [k] at the pole position gives the ap- 
proximate location of the resonant peak in the on-shell 
S'-matrix and the Im [k] gives the width of the peak; the 
Q-value of the resonance is just Q — 2Re [k] /|Im [k] |. 
These statements hold for k z = resonances; for the spi- 
ral modes there are solutions for all k z and for discrete 
complex values of the reduced wavevector 7. Once the 
resonant wavevector has been found the coefficients of 
the outgoing waves can be determined to give the emis- 
sion pattern and polarization properties of the emitted 
radiation. As is well-known, the Hankel functions with 
complex k and large argument grow exponentially and do 
not provide normalizable fields at infinity. This is unim- 
portant for studying the emission patterns as a function 
of far-field angle or polarization properties, although it 
may cause some practical difficulties in numerical algo- 
rithms. If desired, this unphysical feature of the solutions 
can be avoided for a given resonance by adding a tunable 
imaginary part of the index of refraction to yield a so- 
lution with real k outside the dielectric. This imaginary 
part represents linear amplification in the medium and 
would give an estimate for the lasing threshold for that 
mode if mode competition effects were negligible. One 
finds that these real k solutions are continuously related 
to the complex k quasi-bound states at real index and 
have approximately the same spatial properties except 
for the absence of growth at infinity. 
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Fig. 3. A comparison of scattering and emission pictures 
for quasi-bound modes. The complex quasi-bound mode 
frequencies are plotted on the Re [kR] — Im [kR] plane. 
On the back panel we plot the real k S'-matrix, scatter- 
ing cross-section at 170° with respect to the incoming 
wave direction. Notice that the most prominent peaks in 
scattering intensity are found at the values of k where a 
quasi-bound mode frequency is closest to the real-axis. 
These are the long-lived resonances of the cavity. Also 
visible is the contribution of resonances with shorter life- 
times (higher values of Im [kR] ) to broader peaks and the 
scattering background. Calculations are for a dielectric 
cylinder with n — 1.5 and for k z = 0. 



3. Polarization of spiral modes in the far-field 

Having outlined how to solve for the quasi-bound spi- 
ral modes we now analyze their polarization properties. 
Strictly speaking polarization of the time-harmonic elec- 
tromagnetic field cannot be defined inside the cavity or 
in the near field around it as the electric and magnetic 
fields need not be perpendicular to one another or to the 
direction of energy flow. In the far-field on the other 
hand, where the radiation is well-approximated locally 
by a plane wave, we should be able to analyze the po- 
larization of the emission from spiral modes in conven- 
tional terms. Assuming the matching problem is solved, 
the coefficients v ln ,(^ m determining Ef{p, ip), Bf(p, ip) 
are known and these fields can be differentiated in or- 
der to find E± everywhere outside the rod. These re- 
lations simplify if we use the large argument expansion 
of the Hankel functions and their recursion relations to 
find the relative magnitudes of the 2?-field components 
in cylindrical coordinates as p — ► oo , 
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(9 









'k z /j 2 J2 V n 
k/l2 E Cn 



im{ip- 
im{tf>- 
im(ip- 
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fhh 2 E z \ 











(13) 



To extract the polarization at a particular angular direc- 
tion ip we need to recognize that far away from the rod 
the radiation is not propagating in the radial (p) direc- 
tion with respect to the cylindrical coordinates centered 
on the rod axis but is instead propagating at angle a be- 
tween the p and z directions determined by Snell's law 
for the z motion (See Fig. We can then rotate our 




sec aB z | (1.4) 

, sec aE z 



In this rotated coordinate system in the farfield, the elec- 
tric field (on the RHS) only has two components in the 
plane transverse to the propagation direction. Thus, 
the polarization of the electric field in the far-field is 
then determined by the ratio of these two components, 
E z (ip)/B z (ip). If these two field amplitudes have zero 
phase difference we have linear polarization in a certain 
direction which can vary as the angle of observation ip 
varies; if there is a non-zero phase shift A between them, 
then we typically have elliptical polarization except in the 
specific case of A = n/2 and \B Z \ 2 = \E Z \ 2 corresponding 
to circular polarization. It should be noted that when 
the angle a is complex k z > k then the wavevector 72 
of the outgoing Hankel functions is pure imaginary and 
there is no propagating radiation as p — > 00; this cor- 
responds to the total internal reflection condition with 
respect to the z motion of internal ray in the cylinder, 
nsin6* = 1 (see Fig- EJ , for which there is no evanescent 
escape. It should be emphasized however that this does 
not correspond to the true total reflection condition for 
spiraling rays, which comes at smaller 9 except in the 
case of normal incidence in the transverse plane. In par- 
ticular one can ask about the polarization states in the 
far-field of spiral whispering gallery modes, which emit 
solely by evanescent escape. 

The analysis up to this point has been exact for a 
dielectric rod of arbitrary cross-section and the vari- 
ous formulas can be used to solve for both the reso- 
nance wavevectors and polarization properties of the spi- 
ral modes of such a system. We have developed and 
implemented a numerical algorithm to do thisii; we in- 
tend to describe the algorithm and present results for 
deformed cylinders in a subsequent paper^ At this point 
we specialize to the problem of spiral modes of a cylinder 
for which a number of analytic techniques are possible 
which will allow us to develop a useful physical picture. 



4. Quasi-bound resonances in the cylinder 

Focusing now on the case of the cylinder (circular cross- 
section) we have the immediate simplification that the 
Helmholtz equation and boundary conditions separate in 
cylindrical coordinates and we can look for solutions cor- 
responding to a single component of the angular momen- 
tum, m, instead of the sums in Eq. 1111 In the following 
we will use the following notation: 



J m ■= J m (liR) and H m := H+(j 2 R) 



(15) 



where 7, = \Jn 2 k 2 — fcf, i € {1,2} and R is taken on the 
boundary of the domain. With this convention we can 
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write the ansatz for the cylinder 



Here, j is the radial mode index enumerating the solu- 
tions for a given m. Using the boundary conditions for 



J 



the continuity of the field, Eqs. @ and (JSJ), we get the 
relations 



E<(r: 


m,j 


) = am j m (cf^ry m * 


r < R 


(16) 




Bf(r; 


m,j 




r < R 


(17) 




E>(r: 


m,j 




r>R 


(18) 




B>(r: 


m,j 




r>R. 


(19) 





C - ^ULf 

— „ Cjri 



(20) 



Using this, Eqs. © and Q can be rewritten in the fol- 
lowing form 



n 2 cos 2 a H m d p J m 



n 3 ) sin 8 J m H n 



n 2 cos 2 9 J m d p H m 



cos 2 a H m d p J m 
im(n 3 — 



— n 2 cos 2 9 J m d p H m 
n) sin 9J rn H rn 



0, 



(21) 



where the angles are given following the convention in Fig. EJB), C), with tan0 = k z /ji. 9 is the interior angle of the 
ray spiraling up with respect to the (x, y) plane and a the corresponding exterior angle. In order for this system to 
have a non-trivial solution the determinant needs to vanish, resulting in: 



(1 -n 2 ) 2 m 2 sin 2 6< 



1 



l 

J- m H>n 



[cos 2 a H m d p J m - cos 2 9 J m d p H m ] 
[cos 2 a H m d p J m - n 2 cos 2 9 J m dpH m ] 



(22) 



'm J - L m 
QTM qTE 



r 



Where wc have defined 
1 



TE 



G 



G™ = 



l 

J m. Hrn. 



[cos 2 a H m d p J m — n 2 cos 2 9 JmdpHm] 
[cos 2 a H m dpJ m ~ cos 2 9 J m dpHm] ■ 



(23) 

This form is useful since the left-hand side is independent 
of k and vanishes as 9 — > for all m, yielding 

= G TE ■ G™ = [HmdpJm - U 2 JmdpHm] (24) 

x [H m dpJ m ~ JmdpHm]- (25) 

The vanishing of the left bracket describes the reso- 
nance condition for the usual two-dimensional TE modes 
and the vanishing of the right bracket that for the TM 
modeSfiSii 3 " so one recovers the correct limiting behav- 
ior as 9 -— > (k z — * 0). In order to find the resonance 
wavevectors for the spiral modes at 9 ^ one needs to 
find the complex roots of Eq. i12_?l) . An example of such 
solutions is given in Fig. 0] where the roots were found 
by the SLATEC routine dnsqe>i4 The series of resonances 
for a given value of m will be labeled by a second integer 
j which indexes the quantized radial momentum in the 
transverse plane. 

Further insight into the solutions can be obtained by 
noting that for k z = and m^O the TE and TM res- 
onance values differ so that at a typical TM resonance, 
for example, the factor G TE will have its typical order of 



magnitude while the factor G™ vanishes. By continuity 
of these functions with 9 we can expect for small 9 that 
the resonances will have one of the factors G™' TE small 
while the other is not, and that the corresponding reso- 
nance will have a TM or TE character, i.e. the electric 
field will be predominantly in the z-direction or predom- 
inantly in the x — y plane. The TE-like resonances will 
then show a large width (imaginary part) near the Brew- 
ster angle while the TM-like resonances will not. As the 
tilt angle 9 increases the spiral modes become full mix- 
tures of TE and TM modes and it is no longer possible to 
classify them in this manner. In Fig.^J 9 is small enough 
to classify them as TE- and TM-like and the coloring 
represents this classification. 

Once the resonances are found, the polarization in the 
far field can be determined by rewriting Eq. (|21|l using 
the functions G™< TE as 

n 3\ 



m (n — n"j sin f 



n 

m (n 3 



2QTE 

— n ) sin ( 



(26) 

The coefficients a m ,^ m determine the ratio of E z to B z 
in the far-field by 



p = 



E> 



sin 9 



n 2 G TE 



(28) 



We will however defer detailed analysis of the polariza- 
tion properties of spiral resonances of the cylinder until 
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Fig. 5. Exact calculation of the coefficient for the reso- 
nant shift (circle) TM-like, (cross) TE-like. The coeffi- 
cient is expected to be unity for small sinx. 



Fig. 4. TE like resonances (red circles) and TM like res- 
onances (blue circles) for a cylinder with n — 2,9 = 0.2, 
dashed line corresponds to the critical angle defined by 
71 = ran and the solid line to the Brewster condition 
71 = m\Jn 2 + 1. Note that the TE or TM association 
does not work around the Brewster angle. 

we have developed the theory in a more intuitive semi- 
classical approximation in Sections El below. 

A. Small 9 expansion 

In the limit of small k z , or 9, we see that Eq. ()22l) is of 
order 9 2 on the left hand side. Close to a TE or TM 
like resonance at k = ko,k z = 0, one of the factors G TE 
or G™ is small while the other is not. We expand the 
small term (which vanishes on resonance) to lowest order 
in 9 2 and insist that the resonance condition is satisfied 
at a slightly shifted value of the resonance wavevector, 
nk = n(ko + Afco). One can then show that 

Ak 1 9 . . 

ko 2 

where a is given by a ratio of Bessel functions, which is 
plotted in Fig. 03 Note that the coefficient aw 1 for rela- 
tively small sin x; exactly this value follows from the EBK 
quantization formula near normal incidence discussed in 
Sectional This result can be compared to experiments 
done by Andrew Poon, 6 where a tilted optical glass fiber 
was illuminated with an unfocused Gaussian beam. A 
simple "wavefront matching" argument for this blue shift 
was given in Ref. and is reproduced in Fig. El A n inci- 
dent plane wave propagates along the Ai a b direction and 
is incident onto a rod tilted by an angle a. An upward 
propagating spiral resonance requires that the spiral wave 
be in phase with the incident wave farther up the tilted 
fiber. However upper incident wave of the same phase 
front must travel an extra distance d to reach the fiber 
when it is tilted. Thus the phase-matching condition be- 
tween the internal wave and the extended incident wave 
reduces the effective cavity length for the spiral wave by 
a distance d/n, where n is the refractive index. Fixing 



the quantum numbers of the resonances, this implies a 
quadratic blue-shift as a function of the tilt angle 9. 



Zlab 



A) 




2™ 

B) 

Fig. 6. A) Schematic of wave- front-matching argument. 
Internal spiral wave of a tilted optical fiber with respect 
to the incoming wave. The phase-matching condition 
between the spiral mode and the external incident wave 
reduces the effective cavity length for the spiral wave by 
a distance d/n, where n is the refractive index. B) The 
spiral quadratic blue shift can be interpreted by unwrap- 
ping the circular fiber. The dashed lines indicates the 
wavefront. The wavefront matched path is only 2ira cos 9, 
therefore the resonances are quadratic blue shifted with 
the tilt angle. The figures are adapted from Poon£ 

In Fig. [7| we compare the blueshift obtained from the 
exact numerical solution of Eq. \\X2\\ to the small 9 expan- 
sion Eq. I|29|) . The agreement is quite good. We will see 
at the end of the next section that the EBK quantization 
will give the same quadratic blue shift. 
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the vector Helmholtz Eq. © can be written as 

= (r) = A ie ^ Sl ^ + A 2 e i ~ <S2{r) (30) 
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Fig. 7. Shift of the resonance condition with the quantum 
numbers j = 21, m = 20 (sinx « 0.33) and n = 2 with 
respect to the internal tilt angle 6>. Crosses are the nu- 
merical solutions following Eq. (|32|l , and circles following 
the small 9 expansion Eq. H29fl . Note the discontinuity 
in the exact numerical values around 8 = 0.33; this is 
due to the onset of the polarization critical angle which 
will be discussed in Section This discontinuity occurs 
beyond the regime of the small 9 expansion. 



5. The EBK quantization conditions 

As already noted, using the quantized solutions of res- 
onance condition Eq. (|22|l . we can easily calculate the 
polarization in the farfield. However, to get insight into 
the dependence of the polarization on the internal ray 
motion, a more appealing approach is to use the eikonal 
method, and the Einstein-Brillouin-Keller (EBK) type 
formulation of the resonance conditions. The eikonal 
method in general refers to finding approximate solution 
of the wave equation with a specific ansatz which is ex- 
pected to be good in the short-wavelength limit (k — ► oo). 
The EBK method describes how to apply that ansatz to 
boundary value problems, typically assuming Dirichlet 
or Neumann boundary conditions*^ The approach has 
been used for example to find approximate quantization 
formulas for the circular and elliptical billiard systems. 16 
Recently it was generalized to treat the scalar Helmholtz 
equation for a dielectric billiard in two dimensions. 10 
However it was also emphasized that the EBK method 
only works for the small subset of boundary shapes for 
which the ray motion within the boundary is integrable^ 
a point which goes all the way back to Einstein's original 
paper m 1917iii The motion of a ray within an infinite 
cylinder is also integrable (the energy and z-components 
of linear and angular momentum are conserved) and a 
generalization of the EBK method should work in this 
case also. The necessary generalization is to introduce 
the boundary conditions appropriate for the coupled E z 
and B z components of the field; this will lead to a gener- 
alization of the Fresnel coefficients for a plane interface. 

We study the vector Helmholtz Equation © for 
E z (x,y), B z (x,y) in the semiclassical limit fc,7i,2 — ► oo; 
in this limit we expect the solutions to have rapid phase 
variations and relatively slow amplitude variations. The 
generalized EBK ansatz for the quasi-bound solutions of 



where A± t 2 are two-component vectors and S is the 
Eikonal. Note that all the functions are defined on the 
two-dimensional x — y plane. Following Refs liotflril we 
can write the general quantization condition 

7 <j> dq-VS = 271-/; + i= 1, 2. (31) 

Here the quantity VS is the gradient of the phase func- 
tions Si , S2 considered as the two sheets of a double 
valued vector field defined on the cross-section, li are 
integers and Ti refer to topologically irreducible set of 
loops. To avoid confusion in this section we temporarily 
drop the subscript 71 — > 7. Keller showed that in order 
for the EBK solution to be single-valued it is necessary 
that these loop integrals of the phase be quantized and 
that any two topologically incquivalcnt and non-trivial 
loops can be chosen. $i is a total phase shift due to 
caustics and boundary scattering; for the scalar two- 
dimensional Helmholtz equation in a circle these phase 
shifts are known for Dirichlet, Neumann and dielectric 
boundary conditions. For the case of a dielectric circle 
the phase shifts are complex in general, representing ei- 
ther refraction out of the circle or the phase shift at the 
boundary due to total internal reflectionpiS In order to 
get the appropriate ray dynamics for the spiral modes it 
is easily shown that the eikonals Si . Si must be identi- 
cal to those of the two-dimensional circular billiard (i.e. 
k z = 0); the new feature here is an additional eigenvalue 
condition on the amplitude two-vector, which will deter- 
mine the eigenpolarization directions. This will lead to 
a modification of the phase shifts <&i with respect to the 
two-dimensional case. 

Before discussing the latter point we briefly review 
the quantization relations assuming the $i are known. 
Two conventional loops for implementing the quantiza- 
tion conditions are shown in Fig. |SJ The first loop, T\, 
goes just outside the inner turning point of a ray of fixed 
angular momentum; VS points in the direction of the 
ray so this integral just yields the length of the caus- 
tic for this ray. No caustic surfaces are crossed and the 
path does not touch the boundary, so the phase $1 = 0. 
Thus the first loop gives a relation equivalent to angular 
momentum conservation, 



where we have replaced the integer l\ — > m to conform 
to our earlier notation. 

The second loop T 2 gives the quantization condition 
for the reduced wavevector 7 in terms of the path length 
L of the loop (the vector field VS is everywhere parallel 
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to this path) 
1 L 



I' 1 



2 cos x — 2 
2irj + $ 2 , 



(f 



X sin x 



1? 



(33) 



where we have replaced the integer I2 by j. This is ex- 
actly the relation we would find for the two-dimensional 
billiard problem, except that the transverse wavevector 
7 = ^Jn 2 k 2 — k 2 has replaced the full wavevector nk and 
the appropriate phase shift <J> 2 needs to be determined, 
and will differ from the two-dimensional case. 




Fig. 8. A) Path of the first curve iy B) Second path r 2 , 
of length L. 



A. Semi- classical boundary conditions for dielectric rod 

To treat the spiral modes of a dielectric cylinder in 
this approach we need to project the three-dimensional 
boundary conditions corresponding to Snell and Fresnel's 
laws into two dimensions. The relevant angles for this 
projection are all defined in Fig. [21 Since we are in the 
ray optics limit we can regard the "scattering" of the 
eikonal from the boundary as the scattering of a plane 
wave from the tangent plane. Because of our assumption 
of only outgoing waves, it is sufficient to assume only an 
incident, reflected and a transmitted wave (see Fig. |5J. 

1. Generalized Snell's law 

The form of each of the fields is given by 



with 



E{i,r,t} (34) 



The gradient of the Eikonal V5 gives the direction of the 
the ray and is of constant length |V5| = n. We find 

-171 cos x d n S =172 cos a 
«7i sin x dtS l =i^2 sin a 



d n S l =iji cosx 
d t S l = iji sinx 



d t S r 



(35) 
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Fig. 9. Schematics of scattering on the projected plane. 
We will expand the wave solution inside into an incoming 
component and a reflected ^> r . The outside transmit- 
ted component is given by ^ . 



The first set of boundary conditions, the continuity of 
the field across the boundary, Eqs. and (JSJ becomes: 



z \ e hiS' 



i-t\S T 



B 



e~"" . (36) 



Since these equations need to hold everywhere on the 
boundary, the phases need to be equal, thus yielding 



7l 5 l = 7l ^ = 7 2 5 t . 



(37) 



Using the fact that the tangent components are continu- 
ous, we obtain 

72 

71 sin x = 72 sin a sin x = — sin a. (38) 

7i 

This equation can be identified as the projection of 
Snell's Law into the transverse plane. Note that when 
k z = we recover the usual result 



n sin x = sm o~. 



(39) 



We can write the projected Snell's Law in a completely 
geometric fashion noting 71/72 = tan 9/ tan a, 



7! . sin a cosy . . . s 

sm a = — sm x = ; — 7. sm x = J (") n sm X- (40) 

72 cos a sm 



With the function f(6) given by 



/(*) = 



cos 9 



yl — n 2 sin 2 9 



1 - sin^ 9 
. — n 2 sin 2 1 



> 1. 



(41) 



Hence, as is clear geometrically, the projected Snell's law, 
leads to total internal reflection of the projected motion 
before the critical angle sinx = 1/n is reached in the 
plane (this is simply because the actual angle of incidence 
is steeper than the projected angle due to the z-motion); 
the function f(9) which determines the effective critical 
angle is plotted in Fig. ^| 
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Fig. 10. Functional dependence of f(9) for n = 2. Note 
that /(£?) diverges at the critical angle when sin# = 1/n. 



2. Generalized Fresnel Matrices 

The kinematics of the projected ray motion has been de- 
termined above simply from the continuity of the tan- 
gential components of the fields; the transport of ray 
flux across the boundary will now be determined from 
the normal derivative boundary conditions. 



23' 



23 r 



i72S* 



(42) 



where the matrices 23 are derived from the boundary 
conditions Eqs. iJBJl and Q and given by the matrices 

fg(i,r) _ f(n- n 3 ) sin 6 ■ d t cos 2 a ■ d n \ 

n 2 cos 2 a ■ dt (n 3 ~ n) sin 6 ■ d n I ^ ' 



23* 





n 2 cos 2 8 ■ d„ 







0, 



(44) 



Here d n and dt are the normal and tangential derivatives. 
We can now relate the incoming field to the outgoing and 
the reflected using the boundary conditions Eq. and 
Eq. B3> 



(45) 
(46) 



where R and T are the general Fresnel matrices given by 
R = (23 r - »') _1 (23* - 23 l ) (47) 

and similarly 



T = (23* — 23 r ) 1 (23 l - 23 r 



(48) 



In the limit 9 — > =>■ a — > these 2x2 matrices become 
diagonal and take the form 



R = 



T 







n cos x~ cos °" 
n cos x+cos o" 

q COS X~ n CQS & 

COS X + n cos t7 



' 2n cos x 
n cos x+cos cr 







2 cos x 



= R = 



T 



r s 
-r t 



(49) 



cos x+ n cos a 



t s 

t„ 



(50) 



We recognize the diagonal elements as the Fresnel coeffi- 
cients for TM (denoted by subscript s) and TE (subscript 
p) plane waves incident on a dielectric interface. The di- 
agonal nature of the matrices implies that for — these 
polarization states are preserved. When 6^0 the matri- 
ces have off-diagonal elements implying the mixing of po- 
larization states upon reflection (strictly speaking these 
matrices mix E z and B z upon reflection, which will be 
shown to be equivalent to rotating the local polarization) . 
Recall that we are applying the EBK quantization condi- 
tion (|31|l to the two-vector [E z {x, y), B z {x, y)); when we 
integrate around the path B in Fig. [S] which touches the 
boundary, we must impose the correct dielectric bound- 
ary conditions on this two- vector. In general this changes 
the ratio of E z to B z and will lead to a multi-valued so- 
lution as we complete the closed loop (for uniform in- 
dex rods, only boundary scattering leads to the rotation 
of the two- vector) . Therefore in order to have a single- 
valued solution the ratio of E z to B z must be unchanged 
upon reflection, i.e. the two-vector must be an eigen- 
vector a of the reflection matrix, 



Ra = Aa. 



(51) 



We thus see that for the spiral modes of the cylinder 
there are two allowed mixtures of TM and TE polariza- 
tion for each resonance labelled by angular momentum 
m and wavevectors 71 , k z ; the eigenvalues of the R ma- 
trix A = e 171 will give the extra phase shift $2 needed 
to complete the EBK quantization condition in Eq. 
($2 = Tj + ir/2 where the term 7r/2 comes from the caus- 
tic phase shift). As already shown above, at 8 = the 
R matrix is diagonal, conventional TE and TM states 
are eigenvectors and the eigenvalues are just the Fres- 
nel reflection coefficients, r s ,r p . These have the familiar 
property of being purely real and less than unity for inci- 
dent angle \ below total internal reflection, and complex 
numbers of modulus unity for angles above TIR. Hence 
one has pure reflection and refraction with no phase shift 
below TIR, and pure phase shift and no refracted wave 
above TIR. The phase shifts for totally internally re- 
flected TM and TE waves are different and well-known 
functions of x~ This familiar behavior is modified for 
the spiral modes. 

The behaviour of the eigenvalues A = e lr) of the ma- 
trix R for k z ^ is shown as a function of the angle 
sin x = m/j in Fig. ^2 Note that rj is in general complex, 
and we are plotting here its magnitude. The magnitude 
of the eigenvalues are different up to a point between the 
Brewster angle (vertical dashed), and the CA (vertical 
dashed). At this point the the eigenvalues become com- 
plex conjugates of each other. We will call this point, the 
polarization critical angle (PCA) and will explain how it 
is related to the far-field polarization below. Hence spiral 
modes acquire a phase shift upon reflection before they 
reach the critical angle. Not until the CA do A^2 lie 
on the complex unit circle as they should for TIR. Thus 
we have a new phenomenon, a phase shift for a refracted 
wave. Having determined the two possible values of <f>2, 
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Fig. 11. Absolute value (solid red and black lines) and the 
phase divided scaled by 7r (dashed red and black lines) 
of the eigenvalues of R for n = 2 and tan (9 = 0.2 vs. 
sinx- The dotted black horizontal line is the Brewster 
angle and the black dashed line is the critical angle. Red 
indicates the TE-like component and blue the TM. 



we can write the quantization condition for the spiral 
resonances as a transcendental equation 



27W 1 — +2m arcsm I — J — 2ttj+— +w/k+Q+i In \r\, 

/ (52) 
where j, m are integers, A = re 1 ^, and we have used the 
angular momentum quantization condition Eq. Ij32(l . We 
can simplify this result further to get an explicit solu- 
tion for rays near normal incidence in the plane, so that 
sinx — > 0: 



7 



m 1 



f(x,o) 



(53) 



where the contributions of the phase shift and the loss 
at the boundary have been combined in /(x)- We will 
analyze the function /(x> 0) further in the next section 
on polarization properties of the spiral modes. 

From Eq. (|53[l we can derive the blue shift for small 9 
by noting the the right hand side without f(x, 9) is just 
the resonant condition for the circle. From the definition 
of 7 = \J nk? — / 



nk(l 



we can write 



/)/,- = nk Q [ 1 + |<9 2 



1 



(54) 



where nk Q is the resonance condition for the circle 9 = 0. 
Comparing to Eq. I|29l) above we see that this implies that 
the coefficient of the small 9 quadratic blue-shift should 
be a = 1 for small sinx, just as we found in Fig.JSJabove. 
As noted there, this coefficient changes slightly when the 
polarization critical angle is reached; this change is cap- 
tured by the contribution from /(x, 9) we have just ne- 
glected. In Tab. [I]wc compare the resonances found by 
the exact wave matching method (Eq. \22l ) and by the 
EBK method finding good agreement for 9 = 0.1,0.2. 

In general we can always calculate the R matrix via 
Eq. I|47|) . find the eigen-polarization directions, and sub- 
sequently act with T on the internal eigen-polarizations 



Table 1. spiral resonances of the cylinder with n = 2, 9 = 
0.1,0.2. We compare the resonances calculated by the 
solution of Eq. lj2"2)l and the EBK method, finding good 
agreement. Although resonances can no longer be classi- 
fied as TM or TE, classification as TM-like or TE-like is 
based on which factor Q TE >™ j s small at the resonance 
as discussed above 
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to determine the corresponding polarization in the far- 
field. A physically more transparent method to do this 
involves the a reformulation of the problem in terms of 
the actual polarization vector, rather than the two- vector 
(E Z ,B Z ). Below, we will introduce an equivalent ma- 
trix description in terms of the Jones Algebra which we 
will subsequently show to be exactly equivalent to the 
(E Z ,B Z ) description. 

6. Jones formulation of polarization properties 

We introduce the parallel and perpendicular components 
of the electric field, E p and E s in the local coordinate 
system defined by the plane of incidence. The Jones vec- 
tor— which describes the local polarization is 



E 




(55) 



with Eo s and Eq p being the magnitude of the electric 
field and the phases (j) s and <p p . The Jones matrix for 
the reflection and transmission at a dielectric interface is 
give by 



K r 

't v o\ 
t J 







(for reflection), (56) 



(for transmission) (57) 



and would describe any series of reflection in the same 
plane of incidence. However a spiraling ray in the cylin- 
der changes its plane of incidence at each reflection, so 
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Fig. 12. Comparison of the of the eigenvalues of the ro- 
tated Jones matrix, absolute value (solid line red/black) 
and the phase (divided by tt) (dashed line red/black), to 
the eigenvalues of R (red/black circles). Parameters of 
the calculation are tan (9 = 0.2 and n = 2. 



that we need to rotate our coordinate system into the 
new plane of incidence between each reflection before we 
apply the Jones reflection and transmission matrices. (It 
is precisely this rotation of the plane of incidence for a 
spiraling ray which is the physical reason behind the non- 
conservation of TE or TM polarization). If the angle 
between the two planes is given by £ we can do this by 
multiplying the Jones vector by the rotation matrix 9t(£) . 
This angle £ can be determined for a general cylindrical 
symmetry^! Due to the rotational and translational sym- 
metry of the cylinder, the polarization at each angle <f> on 
the cylinder must be the same for any value of z and must 
be described by the same Jones vector when referred to 
the plane of incidence at that point. Therefore the Jones 
vector emerging from a reflection must be the same as 
the incident Jones vector once the coordinate system has 
been rotated into the new plane of incidence. This yields 
an eigenvalue condition for the Jones vectors describing 
the spiral modes of the cylinder, 

'ftl—ftL (58> 

where 2 = 9K(£)J r - Thus the two polarization states 
are the two eigenvectors of 3 and their eigenvalues 
describe the phase shifts and refractive losses at each 
reflection, just as do the eigenvalues of R for the two- 
vector E Z ,B Z studied earlier. Therefore, the two matrices 
3 and R are related by a similarity transformation and 
have the same eigenvalues*^ which is confirmed in Fig. 1121 

7. Polarization critical angle 

Realizing that the eigenvalues of interest are obtained 
from the product a rotation and a diagonal Jones matrix 
with known entries allows us to understand the behavior 
of the these eigenvalues rather simply. The eigenvalues 
can now be written in terms of the rotation angle of the 
plane of incidence, £ and the Fresnel reflection coefficients 
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Fig. 13. A) Absolute value of the two eigenvalues v\ t 2 of 
the rotated Jones matrix (solid red and black). Phase of 
the eigenvalues (divided by tt) is plotted in dashes. The 
Eigenvalues become complex at the point where the two 
curves meet and join. This point lies between the Brew- 
ster angle (dashed vertical black) and the effective critical 
angle (dotted vertical black). Calculated for tan9 = 0.2, 
n = 2. B) (Black) The sine of the polarization critical 
angle PCA at which the eigenvalue of R gets complex. 
(Dotted) Sine of the Brewster Angle, (Dashed) sine of 
the critical angle of total internal reflection. 



v\,2 = - cos£ (r s - r. p ) ± - y cos 2 £ (r s - r p ) + Ar p r s . 

(59) 

Consider the discriminant of the eigenvalues, D = 
cos 2 £ (r s — r p ) 2 + 4r s r p , which determines whether they 
are real or complex. Recall that the TM Fresnel reflection 
coefficient, r s is real and positive for all angles below the 
critical angle passing through unity and becoming com- 
plex and unimodular above the critical angle; whereas 
the TE reflection coefficient, r p is real below the crit- 
ical angle but becomes negative at the Brewster angle 
and passes through negative one before becoming uni- 
modular and complex above the critical angle. It follows 
that D will always be positive and vi t 2 real for angles 
of incidence below the Brewster angle. However for any 
non-zero value of the rotation angle £, D will become 
zero before the critical angle since at the critical angle 
D = 4(1 — cos 2 £) is negative. The value of the incidence 
angle £ when D = is the polarization critical angle 
(PCA) which we have already mentioned above. Since 
it occurs when both r s , r p have absolute value less than 
unity, the eigenvalues also have modulus less than 
unity and we have a phase shift upon reflection while a 
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fraction 1 — | z^i , 2 1 2 of the incidence wave is refracted out. 
This fraction can be calculated and does not correspond 
to either of the usual Fresnel transmission coefficients for 
TM or TE. The behavior of the eigenvalues just described 
is shown in Fig. El A); the behavior of the polarization 
critical angle vs. sinx is shown in Fig. EHB). We see 
that for small 8 the onset of the the phase shift is close 
to the critical angle (CA); as 9 varies the PCA moves 
close to the Brewster angle and then for close to the 
CA, where for any sin% we will have TIR, the PCA re- 
turns to the CA. This analysis allows us a simple un- 
derstanding of why there is a PCA which precedes total 
internal reflection. The TM and TE components of the 
spiral resonances have no relative phase shift at reflection 
until the Brewster angle; at the Brewster angle the TE 
component picks up a t phase shift, which mixes with 
the TM component to give a phase shift between zero 
and 7r. Right at the Brewster angle r p vanishes and the 
local TE component of the resonance is filtered out, but 
as the TE reflectivity picks up above the Brewster angle 
this phase shift appears before total internal reflection 
condition is reached. It is interesting to note that due 
to the finite phase aquired at the PCA, the resonance 
condition Eq. i|53|) is slightly changed. The term PCA 
suggests that at this angle the polarization properties of 
the spiral modes change. We shall see that this is the 
case in the next section. 

8. Polarization properties of spiral modes 
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Fig. 14. A) Phase differences between the E z and B z 
field, for TE and TM like modes. B) Absolute values of 
the ratio of the eigenvector components of R. For normal 
incidence sin% = the modes are clearly either TE or 
TM, right at the PCA they completely mix. Calculations 
done for n = 2, tanf? = 0.2. The solid vertical black line 
is the effective critical angle, the dashed line the effective 
Brewster angle. 



When the eigenvalues of the Z matrix are real, the eigen- 
vectors can be chosen real (up to an overall phase) and 
thus there is no relative phase shift between E s and E p . 
While in general polarization is not well-defined inside 
the dielectric, in the short wavelength limit we are now 
examining it is well defined with respect to the considered 
ray direction, and zero relative phase shift implies linear 
polarization of the resonance with respect to the spiral- 
ing ray direction. Above the PCA, when the eigenvalues 
are complex, the eigenvectors also become complex and 
there is a relative phase shift between E s and E p cor- 
responding to elliptical polarization of the internal field. 
This picture is confirmed by the calculation of the eigen- 
values plotted in Fig. El A), B). At small sinx we have a 
large ratio between the eigenvector components and zero 
phase shift, corresponding to linear polarization with the 
electric field almost completely in the z or transverse di- 
rection (TM-like and TE-like). After the PCA we get 
a phase shift approaching tt/2 and equal ratios leading 
to circular polarization right at the CA; above the CA 
we have in general elliptical polarization for the internal 
fields with a calculable ellipticity and with a phase dif- 
ference of tt/2 as expected for whispering gallery modes. 

The far-field polarization arising from this internal field 
can be found by applying to the relevant eigenvector of 
3 the transmission matrix J t and then rotating the re- 
sulting vector by an angle which projects the Jones 
vector onto the the plane perpendicular to the outgoing 



ray direction discussed in Section [S] above. This angle 
can be determined by straightforward geometric consid- 
erations which we omit here and simply state that the 
the far-field polarization vector (E' s , E' p ) is given by 

£R(-e)Ji|a) (60) 

where 9t is the the rotation just mentioned and |a) is 
either of the eigenpolarization vectors of Z- 

In Fig. ll5l A). B), we compare the far-field polarization 
states for spiral resonances of the cylinder as predicted by 
the exact and geometric optics approach. We find good 
agreement between the methods, although the exact so- 
lutions smooth the abrupt behavior near the PCA pre- 
dicted by the geometric optics approach. Above the CA 
we have only evanecent emission, but the eccentricity of 
the two-vector is still finite. Note that as formulated, the 
Jones approach gives a continuous solution for the two 
polarization states, whereas the exact resonances are dis- 
crete and correspond to discrete allowed angles x which 
can be also found through the EBK approach. The Jones 
approach provides a smooth and fc-independent formula 
for the polarization states which agress with those of the 
resonaces. 

9. Summary and conclusions 

We have reduced the Maxwell's equations for a dielec- 
tric rod of arbitrary cross-section to a vector Helmholtz 
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Fig. 15. Respective A) phase differences and B) abso- 
lute values of the ratio of the eigenvector components of 
R (green and red solid) in the farfield. The black cir- 
cles are the exact numerical solutions following Eq. IL'l'li 
(to € [0,50] and 71 < 50). For n = 2, tan0 = 0.2. The 
solid vertical black line is the effective critical angle, the 
dashed line the effective Brewster angle. 



equation for a two-component vector field living in the 
two-dimensional cross-sectional plane. We have devised a 
formulation of the resonance problem for the quasi-bound 
modes (spiral resonances), which can be implemented nu- 
merically for a general cross-section, and shown how the 
polarization state of the resonances in the farfield can 
be determined. Calculations were reported for the case 
of a circular cross-section (cylinder) and the results were 
compared to ray-optical results from an EBK formula- 
tion of the resonance problem in the semi-classical limit. 
We have analyzed the polarization state of the spiral res- 
onances both inside the cylinder and in the farfield, and 
related its properties to the internal ray motion. It was 
shown that as the tilt angle of the spiraling ray with 
respect to the cross-sectional plane is increased, there 
exists a polarization critical angle at which the polar- 
ization changes from linear to elliptical both internally 
and externally and this occurs before the total internal 
reflection condition, so the effect can be measured read- 
ily in the far-field. The physical picture we developed in 
terms of the Jones polarization vectors was useful in un- 
derstanding the PCA and may be useful in generalizing 
the analysis to arbitrary cross-sections for which the ray 
motion can be chaotic. 
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